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Abstract 

Asymmetric electrical conductance is theoretically demonstrated on the surface of a topological 
insulator (TI) in the limit of infinitesimally small forward and reverse biases between two spin 
selective electrodes. The discontinuous behavior relies on the spin-momentum interlocked nature 
of TI surface electrons together with the resulting imbalance in the coupling coefficients between 
the electrodes and TI surface states. The analysis is based on a transmission matrix model that, in 
combination with a phenomenological treatment for the diffusive limit, accounts for both ballistic 
and scattered paths simultaneously. With the estimated conductance asymmetry over a factor of 
10, implementation in the ratchet-like applications and low-voltage rectification circuits appears 
practicable. 
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At the very dawn of semiconductor electronics, the two-terminal rectiher or diode was 
known as the device that breaks the conductance invariance with respect to the electric 
current (or bias) reversal. The signihcance of this nonlinear property along with the simple 
physical implementation has made major impacts on all areas of electronics including the 
information processing, communication, and power systems. On the other hand, the recti- 
hcation effect in a conventional diode, although well dehned under an applied voltage down 
to tens of mV, attenuates with a smaller signal. The current-voltage relation eventually 
recovers symmetry across the zero bias; i.e., GiV) ^ G{—V) as V —?■ 0, where G denotes 
the conductance at the channel voltage V. 

The physical reason for the symmetric low-voltage conductance in the non-coherent 
regime is based on the principle of detailed equilibrium, i.e., the invariance of scattering 


characteristics with respect to the exchange of incident and scattered particles [l| . In quan¬ 
tum mechanics, this property is attributed to the Born approximation that is often treated 
as a universal quality of electronic scattering j^. More formally, it follows the Onsager re¬ 
ciprocal relation G(V, H) = G{—V, —H), where H is the magnetic held. The involvement of 
magnetic held clearly opens a door to the possibility of breaking the conductance symmetry 
with respect to the voltage alone as G(V, H) does not have to be equal to G(—V, H). Sim¬ 
ilarly, the correlation between the magnetic inhuence and the electrical transport may, in 
principle, lead to the asymmetry even for V —)■ 0. Nevertheless, accessing the zero-voltage 
inequality remains a signihcant challenge despite numerous attempts js-JJ]. A system with 
high intrinsic asymmetry could be far more advantageous in realizing the desired character¬ 
istics. 

The spin-orbit interaction (SOI) is thus a primary consideration. Combined with spin 
polarized injection, it was shown theoretically that the SOI can cause transmission im¬ 
balance (zero-bias) in an asymmetric open structure by diherentiating the passes for the 
spin-up and spin-down electrons in the coherent quantum mechanical regime 13|]. Even 
under more conventional transport conditions, spin injection at an arbitrary small voltage 
creates a non-equilibrium spin distribution in the non-magnetic medium. The SOI trans¬ 
fers spin polarization to the electronic flux, directing it according to the spin and not to 
the voltage polarity [l^. This net effect of spin-dependent flux on the total current could 
lead to the zero-voltage conductance discontinuity if the spin polarization is conserved. A 
prominent example of the spin-current conversion is the photo-voltaic effects that appears 
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due to the spin polarization of the photo-excited electrons (e.g., see Ref. llj and refer¬ 
ences therein). However, the directional selectivity of conductance (and its zero-voltage 
discontinuity) of spin-polarized electrons has not been pursued vigorously in conventional 
semiconductors as the effect of scattering asymmetry is often too small to overcome the 
background spin-independent events. Moreover, the spin-related asymmetry decays quickly 
once the thermal energy becomes dominant over the spin-orbital splitting induced by the 
Dresselhaus or Rashba effects. In any case, the intrinsic imbalance relying on the spin de¬ 
pendent properties does exist, making systems with large spin-orbital coupling potentially 
far more advantageous in realizing the desired characteristics. 

In this regard, the recently discovered 3D topological insulators (TIs) provide a uniquely 


promising candidate 


15|. Unlike the traditional semiconducting 2D channels, the electron 


spin on a TI surface is inherently locked to its momentum. Thus, the electron spin po¬ 


larization follows the surface current and vice versa 
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18| . The subsequent directional 


selectivity should immediately result in the ratchet effect when combined with the ferro¬ 
magnetic (FM) contacts that enable spin polarized injection and collection. Indeed, a recent 
measurement illustrated a difference in the TI surface resistance between two FM electrodes 
at hnite positive/negative voltages 1^. Yet, beside being small, it is highly unclear whether 
the observed inequality persists at zero bias. In the ballistic regime, the transition matrix 
symmetry based on charge conservation imposes the conductance reciprocity in a strictly 
two-terminal structure 20|]. At the same time, theoretical calculations in a long channel 
(i.e., the diffusive transport) also reveal the reversible conductance as the initial net po- 

n n 

larization quickly undergoes relaxation after a number of scattering events |21l . I22| . The 
desired asymmetric response clearly requires an additional consideration while not causing 
a signihcant loss of electron spin polarization. 

In the present work, we propose a solution to the long running challenge by exploiting 
the advantages of a TI based structure for highly asymmetric conduction at arbitrarily 
small forward and reverse biases. The key variation from the conventional two-terminal 
assumption, as shown in Fig. [H is the hnite size of the FM electrodes in the driving direction 
(i.e., x), which are also sufhciently away from the outer edges of the TI sample (to avoid 
the boundary scattering). This enables the electrons on the TI surface to establish a non¬ 
isotropic, non-equilibrium distribution following the spin injection conditions of the FM 
contact (with magnetization M along the y axis) that remains signihcant over the scale 
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of spin relaxation length at an infinitesimal channel voltage. In a sense, the TI layer may 
be considered as a spin sensitive scatterer that separates the electrons in the real space 
according to their spin momenta—analogous to that in the Stern-Gerlach experiment. Thus, 
spin specific electron transport can be manifested in the directional conduction so long as the 
polarization survives the relaxation in the channel, leading potentially to GiV —?■ 0“, M) 7^ 
G{y —)■ O’*", M). When the contact width is much larger than the spin relaxation length, the 
electron distribution on average loses the asymmetry and, hence, the polarization specific 
response. In fact, adoption of narrow FM contacts and subsequent electron transport in the 


spin-momentum interlocked channel is rough 
with spin polarized injection described above 


y analogous to the asymmetric open system 
13j, even though the operating conditions are 


drastically different (i.e., non-coherent vs. coherent). 

The origin of the asymmetry between the opposite bias polarities is more clearly illus¬ 
trated in Figs. [2])a) and[2](b). Here, the light (orange) and dark (purple) arrows denote the 
ballistic and scattered paths, respectively. In the low resistive case [Fig. [2][a); forward bias], 
the electrons injected from the FM contact on average possess the TI surface momenta that 
are naturally aligned with the applied electrostatic force. Hence, the conduction is primarily 
via one spin channel directly connecting two electrodes (for both ballistic and non-ballistic 
transport). Once the sign of the applied voltage is flipped [Fig. [2](b); reverse bias], however, 
the contribution of the direct channel involves the minority spin injection that has a signif¬ 
icantly limited capacity and the paths involving spin/momentum altering transition must 
be accounted for. The nature of this imbalance is expected to persist when the bias tends 
asymptotically to zero [see, for instance. Fig. [2][c) as well as the related description for nu¬ 
merical estimation given later in the discussion]. Nevertheless the asymmetric conductance 
does not yield a non-zero current at zero voltage. 

Our theoretical analysis predicts the presence of inequality even at infinitesimally small 
biases, leading to potential discontinuance in the channel conductance as the separation 
between the two electrodes (L) shrinks. In the experimental implementation, the required 
dimension may not be as demanding as it first looks. While the measurements of A on 
the TI surface suggest only about a few nm at the moment, this is due primarily to the 
poor sample quality. Rather, the projected values limited by the intrinsic electron- pho non 
interaction could reach several pm once unintended imperfections can be minimized 23l-l25[. 


Hence, achieving a high-quality material with the critical length around 100 nm appears 
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entirely feasible. 

Before examining the details, it is important to note that the hnite sizes of the contacts 
and the distance between them complicate accurate evaluation of the electrical conductance. 
More specihcally, transmission via both ballistic (direct) and scattered (sequential) paths 
must be accounted for simultaneously even when L < X. As described earlier [see Fig.[2][a,b)], 
this is because some lucky electrons injected with a ’’wrong” spin (thus, a momentum in the 
’’wrong” direction) can still turn around and complete the conduction despite the restriction 
on the 180° back scattering. Hence, the widely used treatment based on the Landauer- 
Biittiker approach is not adequate, requiring further development as discussed below. 

The theoretical formulation starts with the description of electrons in the FM electrodes 
and the TI surface channel. In the actual realization, an atomically thin tunnel barrier may 
separate them (i.e., the TI and the electrodes) for highly efficient spin polarized injection 
as well as to avoid the modihcation of the TI band structure. By considering spin-up and 
spin-down components explicitly, the eigenfunctions in the electrodes can be written as: 


x/TT^ 



( 1 ) 


where s symbolizes the spin index and M(r) is the spatially varying function that is de¬ 
termined by the structure geometry and the constituent material properties. Accordingly, 
= s, where s = ±1 corresponds to M || ±y, respectively. On the TI surface, 
a given spin polarization determines the electron momentum due to the spin-momentum 
interlock. The two spinor wavefunction on the TI surface is thus described by its wavevector 

k = {k^, ky) as: 



w{r) 


hvpiifikx + ky)/E 
1 


exp[i{(3k^x -F kyy)], 


( 2 ) 


where w{r) is the Bloch function (that decays exponentially from the surface along the z 
axis) and vp is the Fermi velocity of TI surface electrons (= 4.5 x 10® m/s for Bi 2 Se 3 ). When 
the electrodes are elongated along the y axis as shown, the invariance in this direction can be 
assumed, enabling the states to be identihed by the energy E and transverse wavevector ky. 
As such, kx simply becomes \J(E/hvp)'^ — ky with the factor (3 denoting the -|- and — signs 
for the rightward and leftward traveling components, respectively. Hence, the TI electrons 
with a given energy and transverse wavevector can be expressed as a linear combination 
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of the spinor eignefunctions. While ky can take both positive and negative values, the two 
components (±) of the x directional momentum are considered explicitly. 

Then, electron transmission between the FM contact and the TI surface states can be 
treated via the tunneling process as in the scanning tunneling microscope 26|. For instance, 
the stimulated surface states astride the narrow FM electrode at a; = 0 can be described as 
Tti = where 9{x) is the Heaviside step function and the coefficient rep¬ 
resenting electron injection in each mode simply becomes ~ AV;, 

if a scattering potential AVJ, is adopted as a formal representation for the tunneling Hamil¬ 
tonian Ht through a uniform barrier. A straightforward calculation subsequently leads to 

E — shvp^fikx -|- iky) 


a^(k) = To- 


2E 


( 3 ) 


where tq results from AVj, and the overlap integral of the spinor-independent parts, u{r) and 
t(;(r). The explicit dependence of Eq. (3) on s and jS clearly illustrates the influence of spin- 
momentum selection; i.e., the imbalance between two transport directions when the contact 
is spin specihc. These complex terms originating from the spin-orbit interaction make the 
expression non-analytic (in terms of Cauchy-Riemann conditions in the 2D parameter space) 
that, if persists in the TI channel, can give rise to the discontinuity in the zero-voltage 
conductance. As an example, the normalized directional dependence of |a^(k)p is plotted 
in Fig. [2](d) for the case of s = —1 (i.e., with the —y magnetization). Evidently, the injected 
distribution favors the initial momentum along the +x direction with approximately the 
80:20 ratio between the majority and minority populations. This preference is independent 
of the electron energy E. The absorption of TI electrons by the collecting FM electrode 
can be handled similarly. The only difference is a slight modihcation in the Heaviside step 
function for Tti as the (3 = + (—) mode now approaches the contact from the left (right). 

Following the overlap integral, the direct conduction through the ballistic paths be¬ 
tween the two electrodes can be evaluated by treating the TI surface state as the inter¬ 
mediary that couples the FM states. This results in the transmission coefficient ttiE) ^ 

where the proportionality to the energy originates from the density 
of modes on the TI surface 27(]. While numerical evaluation of the transmitted flux (i.e., 
|t^P) is difficult due to the uncertainty in tq, it is outright obvious that the relative ratio 
of the conductance can be obtained as the polarity of the magnetization and/or the driving 
bias change. Assuming both electrode are magnetized along the —y axis (s = —1), for in- 
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stance, a switch from the forward to reverse bias (/9 = H —> —) yields conductance variation 
as large as 70 to 1 [more precisely, (1 + vr/4)^ : (1 — vr/d)^]. For convenience, this factor 
illustrating the dependence on the magnetization and bias polarities is defined as 

M _ (1 -g/^7r/4)^ 

“ 2(l + 7rVl6)’ ^ 

The denominator normalizes the directional contribution to the sum of both spin states that 
corresponds to the case of nonmagnetic electrodes. Hence, the actual flux can be obtained 
formally from the product of and the spin-independent, ballistic path conductance (de¬ 
noted as Following the schematic given in Fig. [H we set s = — 1 in the ensuing 

discussion. The case for the +y magnetization (i.e., s = -|-1) can be evaluated similarly. 

Along with the ballistic conduction between two FM electrodes, it is necessary to account 
for the role of the sequential channel involving the scattered paths even when L < A. Due 
to the finite width of the contact region, the electrons initially injected with the momenta 
moving away from the collecting contact can turn around after a few scattering events and 
provide a non-zero contribution. It is particularly relevant if the conduction is in the high 
resistive (i.e., reverse biased) state as indicated in Fig. [2][b). When L is large, on the other 
hand, the non-ballistic current becomes dominant and the channel conductance converges 
to the spin-independent diffusive limit. Unfortunately, a rigorous treatment of this complex 
transport problem is very difficult to attain. Instead, we approach it by estimating the 
effective length of the non-ballistic paths where the scattering events are modeled as elastic 
processes. Typically, the additional distance that an electron injected into the ’’wrong” spin 
state (thus, the momentum) travels away before experiencing the turn around is of the order 
of the mean free path A. Given that the back scattering on the TI surface is suppressed, the 
average detour qX is somewhat longer than A with a value around g ~ 4 following a simple 
analysis. A similar process is also encountered in the absorbing electrode, causing additional 
differences in the paths. Of the two conduction conditions, the high resistive state is again 
expected to experience stronger reflection due to the spin mismatch and, thus, a longer 
path on the collecting end [see Fig. [2][b)]. We neglect this extra spin (or polarity) selective 
contribution as it is not crucial in the main focus of the investigation (i.e., demonstrating 
the nonreciprocal conduction). The resulting approximation tends to underestimate the 
asymmetry ratio; hence, a rather conservative assumption. 

Then, the non-ballistic surface channel conductance can be given as shown on the right 
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side of the expression [2 


Gxi^ti — G^i X 


P+(1-P)(TT^> 


forward 


( 5 ) 


(1-P)+P(T+^’ 


reverse 


Here, G^i represents the (bulk) TI surface conductance of length L (hence, the conductivity 
of LG^i). In addition, the factors p (~ 0.8) and 1—p account for the uneven distribution of 


the injected electrons between the +kx and —k^ components as discussed earlier in Fig. [2](d). 
The portion grouped in the curly parenthesis is subsequently dehned as Tpi, highlighting the 


normalized spin/polarity-dependent component. To obtain the contact-to-contact transport 
characteristics, the (sequential) tunnel resistance between the electrode and the TI (1/G^j) 
must also be considered. Hence, the total resistance in this case becomes 1 /Gti^ti + 2/Gti- 


In a steady state with a nonzero bias V, the total current is obtained by summing the 
contributions from both the ballistic and non-ballistic transport. Thus, the resulting total 
conductance [= I{y)/V] follows the form: 



( 6 ) 


As shown, the exponential decay is considered explicitly in the ballistic contribution (the hrst 
term), while the second accounts for the diffusive conduction described directly above. Under 
a small voltage, the I-V curve is expected to be linear in each of the forward and reverse 
bias regions but with different slopes [i.e., Gtot(U) ^ Gtot(—U); see also the illustration in 
Fig. 2(c)]. Such a possibility is evident from the dependence of Eqs. (jl]) and (jS]) on the 
relevant parameters. 

For the quantitative results, a parametric study is carried out as the numerical values of 
^Ti established. To simplify the analysis, the surface conductance Gti 

dehned for a hxed length A (electron mean free path) is used hereinafter; i.e., ^ f ^ti 

in Eqs. (jS]) and ([6]) showing explicitly the dependence on L. In this regard, the channel 
dimension is described in units of A as well. Figure[3](a) shows the conductance ratios between 
the forward and reverse biases of inhnitesimal amplitude for three different combinations of 
and G^/G^^i. The rapid decay with the increasing L conforms to the relaxation of 
the spin selective injection through the interaction with the environment. Once L S> A, the 
transport properties converge to the drift-diffusion equations for all cases under consideration 


2l| . In the other extreme, the ballistic conduction dominates and a large asymmetry over the 








factor of 10 can be clearly achieved as desired. For the strong rachet effect, the calculation 
results also indicate the need for a large and/or a small G^j in reference to Gti- 

The impact of G^j/Gti and GJ^/Gti is analyzed in greater detail with L hxed at 3A in 
Fig. m^b). The star on the upper left corner indicates the limit of approx. 70 set by T^. 
Evidently, the ratio approaches to the maximum when the conductance is dominated by the 
ballistic transport (i.e., a large GJ^/Gti). As for the dependence on G^i/Gti, the asymmetry 
ratio is not a monotonously decreasing function despite the initial impression. At the low 
limit, the large tunnel resistance between the electrode and the TI effectively blocks the 
sequential conduction involving the scattered paths, leaving the ratio to follow that of the 
direct conduction via ballistic paths. This diminishes any contribution from Tti [Eq. (E])] for 
the asymmetry. Once G^j increases, the scattered paths start to draw the current and the 
ratio decreases rapidly as expected (i.e., a more signihcant role of the diffusive regime with 
a smaller ratio). When G^j becomes even larger toward the high limit, the asymmetry ratio 
deviates somewhat from the simple decay with the value determined by the interplay between 
GtiTti. One potential complication is that G^j and Gj(j may not be controlled 
independently as both of them are related to the tunneling characteristics. Nonetheless, 
a thinner tunnel barrier is preferred since ensuring a large contribution from the ballistic 
paths is more crucial for high asymmetry. On a side note, the white dot in Fig. 3(b) denotes 
the conditions under which the I-V characteristics in Fig. [2](c) are obtained. The piecewise 
linear relation with a sizable difference in the forward/reverse slopes (a factor of 10) provides 
distinctively unique features compared to the conventional diode with a similar current scale. 

In summary, the asymmetric conductance between two narrow magnetic electrodes on a 
TI surface is predicated. This phenomenon requires that the distance between electrodes be 
comparable to the electron mean free path and the electrodes far from the TI boundaries to 
minimize the relaxation. The spin selective effect of the magnetic electrodes coupled with 
the intrinsic spin-momentum interlock causes the conductance imbalance which is far more 
signihcant than that based on the structural asymmetry. Our theoretical estimate indicates 
a clearly measurable conductance imbalance between the forward and reverse biases even at 
inhnitesimally small voltages. A key for experimental verihcation at room temperature is 
the improved sample quality as the characteristic dimension is limited by A. 

This work was supported, in part, by FAME (one of six centers of STARnet, a SRC 
program sponsored by MARCO and DARPA) and US Army Research office. 
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FIG. 1. (Color Online) A diode-like TI structure with two magnetic electrodes on the top surface. 
Thin tunnel barriers separate the TI surface and the metallic electrodes (not shown). Two red 
arrows denote the magnetization direction of the magnets (FM), while V stands for the applied 
bias. The purple arrows (thicker vs. thinner) indicate the preferred direction of electron flow 
corresponding to the specified magnetization profile. 
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FIG. 2. (Color Online) (a,b) Schematic illustration of electron flow under the forward (low re¬ 
sistance) and reverse (high resistance) biases, respectively. The total current is the sum of the 
contributions from the the direct (orange arrows; ballistic paths) and the sequential channels (pur¬ 
ple arrows; scattered paths). The magnetization of electrodes is assumed to be along the —y axis, 
(c) Anticipated asymmetric I-V characteristics near the zero bias with the rectification ratio of 10 
(red). Also shown for comparison are those of a conventional space-charge diode with a comparable 
current level (black). A related description can be found in the caption for Fig. 3(b) as well, (d) 
Normalized tunneling strength |a^(k)p between the FM electrode (with the —y magnetization) 
and the spin-momentum interlocked TI surface states. The area under the curve is set to 1. The 
result clearly shows the imbalance between the right- and left-moving electron distribution. 
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FIG. 3. (Color Online) (a) Conductance ratio between the forward and reverse biases of infinitesi¬ 
mal amplitude (i.e., the conductance asymmetry or discontinuity) as a function of the normalized 
channel length L/X for three different combinations of (G^j/Gti, GJ^J/Gti). (b) Contour plot of 
the conductance asymmetry ratio in the G'ti-G^/G^ i parameter space for L = 3A. The star 
indicates the limit of approx. 70 that is set by the purely ballistic transport. The numbers denote 
the calculated asymmetry values. The lighter (darker) color signifies the higher (lower) ratio. The 
white dot indicates the parameters used for the I-V calculation (spin diode) shown in Fig. [2)(c). 
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I. CALCULATION OF BALLISTIC TRANSMISSION COEFFICIENT BETWEEN ELECTRODES 


With the only assumption of spin conservation, the coupling between the electrodes and the TI surface [i.e., 
{tp^^{E,s)\AVb\'tpJ^{kx,ky))] follows the spin selective rule, which results in the coupling coefficient dependent on 
the TI electron momentum as shown by Eq. (3) in the main paper. Here, it should be noted that the states in the 
electrode is described by the energy E and the spin polarization s, while that in the TI is by the momentum according 
to the linear dispersion relation and the spin-momentum interlock. The overall physical picture for the direct ballistic 
coupling between two electrodes can be constructed by treating the TI surface states as the intermediary. Then, by 
taking advantage of the overlap integral through the tunnel barrier mentioned above, a formal expression for this 
coupling coefficient between electrodes can be written as: 


t^iE) = 


E E 

k' K, 


E-E' + i0 




( 1 ) 


where the subscripts for the FM states represent the electrode on the left (L) or the right (K) and (/3fc^,fcy) is the 
momentum of the TI surface state that couples to the electrodes. As in the main paper, we set k'^ a positive value and 
use P to denote the direction of electron flow. Given a width W for the channel, ky is then quantized as ky = rmrjW, 

where m = 0, ±1, ±2, • • •. If the electron has an energy E' = hvF-\Jk'J^ + k'y^, the upper bound for \m\ is determined 
accordingly by E'-, i.e., \m\ < E'W/nhvp- When the width is sufficiently large such that the quantization step tt/W 
is negligible compared to E'/hvp, the summation over all possible ky can be substituted by an integral. The same 
applies to the x direction as well. Indeed, we assume that the x direction is unconstrained, so essentially becomes 
continuous. In addition, it is more suitable to average the result over the system dimension. We adopt the change of 
notation to let denote such an averaged value hereinafter, i.e.. 


(E) = f 

J —( 




E-E’ +i0 


-.dk'dk' 


k'J 




-/2 {,pl^{E)\AVi,\^JHE',e)){i,f{E',0)\AVt\P;l^{E)) E' . 

- - - - -cos 9 -— dOdE . 

_7r/2 E — E' + ifi hvp 


( 2 ) 


This is a quasi-ID treatment that is valid given the invariance over y. The second expression indicates the change 
of variables according to dk',^dk'y = -E—dE'd9 (i.e., the E-k relation). In addition, the phase factor is dropped 
from the equation as it is for the interference of electrons with different wavevectors. The interference phenomenon 
requires very stringent conditions that are generally not supported in the realistic systems at room temperature. The 

k' 

cos 9 factor, or the ^ term, explicitly accounts for the projection along the channel direction. The integral 

over the energy can be carried out by using Gauchy’s residue theory. This is intuitively expected as the resonant case 
with E' = E dominates the contribution. Thus, Eq. m becomes: 


(E) = -iF (V;™(£;)|AW|^J/E,0))(V;Ji(£;,0)|AW|V;™(E))cos0d0. 

d-7r/2 


(3) 


Gombined with the overlap integral given in Eq. (3) of the main paper, the coupling between the two electrodes 
through the TI surface states is finally given as (with = E cos9/hvp and ky = E s\\i9/hvp)' 


{E) = —ITT J 


'^1'^ E — shvpiPk!^ + iky) 


-7r/2 


ko 


2E 


' cos9d9 = —*7r|To/-^ ~ t 

hVF ^ 4 


( 4 ) 
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II. ESTIMATION OF DIFFUSIVE CHANNEL CONDUCTANCE 

As described in the main text, the non-ballistic diffusive conduction can come from the electrons injected with both 
the ’’correct” spin (moving toward the collecting contact) and the ’’wrong” spin (moving away from the collecting 
contact). The electrons in the first group (’’correct” spin) travel the distance of L (i.e., the separation between the two 
contacts) to reach the collecting electrode, while those in the second group (’’wrong” spin) experience longer paths 
(L + q\) since they have to first turn around via scattering events (A being the electron mean free path). When the 
conductivity of the channel is given by LGf^i in the diffusive regime (as specified in the main text), the conductance 
of the former and the latter paths becomes and Gi^iL/{L + gA), respectively, reflecting the difference in the path 
lengths. Then, by treating the probability of an electron injected with the ’’correct” spin to be p' and that with the 
’’wrong” spin 1 — p', the diffusive channel conductance can be written approximately as: 

Gdiff = p'G^i + (1 — —tG^i • (5) 

L + q\ 

When the structure is forward biased, p' corresponds to p (« 0.8) in the main text. For the reverse biased case, it is 
1 — p that matches with p' instead. This leads to Eq. (5) of the main text. Actually, the latter expression adopts an 
additional parameter Tti that is normalized by G^j to highlight the spin/polarity dependence (i.e., Tti = Gdis/G^j). 
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